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ABSTRACT 

Context. NRAO 150 is one of the brightest radio and mm AGN sources on the northern sky. It has been revealed as an interesting 
source where to study extreme relativistic jet phenomena. However, its cosmological distance has not been reported so far, because of 
its optical faintness produced by strong Galactic extinction. 

Aims. Aiming at measuring the redshift of NRAO 150, and hence to start making possible quantitative studies from the source. 
Methods. We have conducted spectroscopic and photometric observations of the source in the near-IR, as well as in the optical. 
Results. All such observations have been successful in detecting the source. The near-IR spectroscopic observations reveal strong 
Ho- and H/3 emission lines from which the cosmological redshift of NRAO 150 (z = 1.517 + 0.002) has been determined for the 
first time. We classify the source as a flat-spectrum radio-loud quasar, for which we estimate a large super-massive black-hole mass 
~ 5 X IO'Mq. After extinction correction, the new near-IR and optical data have revealed a high-luminosity continuum-emission 
excess in the optical (peaking at ~ 2000 A, rest frame) that we attribute to thermal emission from the accretion disk for which we 
estimate a high accretion rate, ~ 30 % of the Eddington limit. 

Conclusions. Comparison of these source properties, and its broad-band spectral-energy distribution, with those of Fermi blazars 
allow us to predict that NRAO 150 is among the most powerful blazars, and hence a high luminosity -although not detected yet- 
y-ray emitter. 

Key words, galaxies: active - galaxies: jets - galaxies: quasars: general - galaxies: individual: NRAO 150 - infrared: galaxies - 
techniques: spectroscopic 



wave and the mm-wave jet, which is, together with the one sid- 
edness of the jet, a clear sign of jet orientation close to the line 
of sight. More intriguing is the evidence of fast "jet wobbling" at 
~1 l°/yr in the plane of th e sky; the fastest rep orted for an AGN 
so far The observations bv lAgudo et alj (l2007h together with the 
cosmological redshift measurement presented in this paper and 
non-contemporaneous X-ray data have allowed to report the first 
quantitative estimates of the basic physical properties of the in- 
ner jet in NRAO 150; i.e., Doppler factor 5 ~ 6, the bulk Lotentz 
factor y ~ 4, the angle subtended between the jet and the line 
of sight (p ~ 8°, and the magnetic field intensity of the flow 
B X 0.7 G. Such large B estimate seems to be compatible with 
the highly non-balistic superluminal motion of the inner jet in 
the source revealed by the mm-VLBI images (fi^pp ~ 3 times the 
speed of light), which have also shown NRAO 150 as a pr ime tar- 
get to study the origin of the jet wobbling phenomenon (lAgudd 

N RAO 150 was first d etected in December 1981 in the opti- 
cal bv lLandau et al.l(ll983h . However, no optical classification or 
distance determination has been reported so far, perhaps due to 
the difficulties to observe the source in the visible range, where 
NRAO 150 is strongly affected by Galactic extinction (Galactic 
latitude b » -1.6°). This problem is partially overcome in the 
near-IR range, where spectroscopic observations from strong 
lined objects can be performed from Earth. Independently of the 



1. Introduction 

NRAO 150j first catalogued by 

iPaulinv-Toth. Wade & Heeschen I (Il966h . is nowadays 
one of the strongest radio an d mm AGN sources 
in the northern sky (e.g. Terasranta et al.l l2005t 



lAgudo. Thum. Wiesemever & KrichbaumI |2010|) . The source 
has be en monitored at cm and mm w avelengths for decades 
(e.g. lAUer et all Il985t iReuter et al.l |I997; Ter asranta etaP 
|2004 and references therein), and has displayed flux densities 
in the range [2, 16] Jy at 2cnu and [1.5,9.5]Jy at 3mrrQ, with 
absolute maxima at the beginning of year 2009. 

At radio wavelengths, on VLBI scales, NRAO 150 dis- 
plays a compact core plus a one-sided jet extending up to 
r i 8 mas with a jet stru ctural position angle (PA) of ~ 30° 
(e.g., iFev & Charlolll2000l) . The fir st set of ultra-high-r esolution 
mm-VLBI images of NRAO 150 dAgudo et al.ll2007h have al- 
lowed to report a large misalignment (> 100°) between the cm- 
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Galactic absorption along the line of sight of the source, the near- 
IR range is the adequate spectral range to detect the strong Ho- 
line f rom AGN at cosmolo gical redshifts between 1.2 and 3.6 
(e.g., iBabbedge et aDl2004 . 

Here, we present the results of our spectroscopic near-IR ob- 
servations, which were successful on detecting, for the first time, 
emission lines from NRAO 150. In Sect.|2]such observations and 
their data reduction procedures are outlined, together with a set 
of optical and near-IR photometric observations performed in 
the 2005-2007 time span. The cosmological redshift determina- 
tion, the classification, the first estimate of the mass of the super- 
massive compact object in NRAO 150 and its accretion rate, as 
well as its broad band spectral energy distribution, are presented 
and discussed in Sect. [3] and H] whereas a summary of our main 
results and our conclusions are provided in Sect. |5] 

2. Observations and Data Reduction 

2.1. Near Infrared spectroscopy 

We have obtained near-IR spectroscopy of NRAO 150 us- 
ing the Long-slit Intermediate Resolution Infrared Spectrograph 
(LIRIS, built at the Instituto de Astroffsica de Canarias 
iManchado etaP 12004 lAcosta-Pulido et al.ll2003h on the 4.2 m 
William Herschel Telescope (WHT) in two epochs (2005 March 
and 2007 January). 

The spectra were obtained using grisms LR_ZJ (bands Z+J) 
and LR_HK (bands H+K). For the observations in 2005, the slit 
width was selected to be 0'.'75, providing resolutions of 700 and 
600 in the Z+J and H+K bands, respectively. For those in 2007 
the slit width was 1" to match the seeing, providing a resolu- 
tion of 500 in the Z+J bands. The slit orientation (PA = 1 18.8°) 
was chosen to include in the slit aperture both NRAO 150 and 
the object 2MASS J03592889+5057547. The observations were 
performed following an ABBA telescope nodding pattern, with 
each AB cycle repeated 3 times. The exposure time for a single 
frame was 600 s, giving a total of 3600 s for both spectral ranges. 
The data were reduced following standard recipes for near-IR 
spectroscopy, using the dedicated software lirisdr, developed 
within the IRAF0 enviroment by the LIRIS team. The basic re- 
duction steps comprise sky subtraction, flat-fielding, wavelength 
calibration, and finally, the combination of individual spectra 
by the common shift-and-add technique . For a more detailed 
description o f the reduction process, see iRamos Almeida et al.l 
(l2006l l2009h . The flux calibration and telluric absorption cor- 
rection are usually determined from observations of nearby stars 
of spectral types AOV or G2V, which are obtained inmediately 
after or before the science observations. In the first epoch the 
selected correction star resulted to be a multiple stellar system 
whose 2MASS near-IR colors indicate a spectral type different 
than AO. We have determined a spectral response function using 
the spectrum of the star 2MASS J03592889+5057547, which 
was observed simultaneously with NRAO 150. The stellar spec- 
tral shape was corrected by a blackbody of temperature corre- 
sponding to an M5 spectral type. A posteriori we derived the 
telluric correction from the normalized spectrum of star 2MASS 
J03592889+5057547. During the second epoch an AOV star was 
successfully ob served. A modifie d version of the Xtellcor rou- 
tine was used dVacca etal.ll2003h to obtain the calibrated flux 
and telluric corrected NRAO 150 spectrum. Relative light losses 

^ IRAF (Image Reduction and Analysis Facility) is distributed by 
the National Optical Astronomy Observatories, which are operated by 
AURA, Inc., under cooperative agreement with the National Science 
Foundation. 



SPECTROSCOPY 



Date 


Spec Range 


Dispersion 


T 

^ e.\p 


Seeing 




yum 


A/pix 


S 


arcsec 


2005 March 23 


0.89-1.33 


6.1 


3600 


an 


2005 March 22 


1.45-2.40 


9.6 


3600 


a.n 


2007 Jan 2 


0.89-1.33 


6.1 


3600 


V'.l 


2007 Jan 4 


0.89-1.33 


6.1 


2400 


1'.'2 


PHOTOMETRY 


Date 


Instrum/Telesc 


Filter 


T»vn (S) 


Seeing 


2005 March 22 


LIRIS/WHT 


J 


400 


1'.'4 






Kj 


400 


1'.'4 


2005 Nov 2 


IAC-80 


V 


9000 


I'.'S 






R 


6600 


1'.'3 






I 


6000 


1'.'3 


2005 Nov 1 1 


IAC-80 


R 


3600 


1'.'4 






I 


4800 


1'.'4 


2006 Feb 03 


Liv. Tel. 


V 


600 


1'.'4 






r' 


300 


1'.'2 






i' 


200 


1'.'2 


2006 Feb 17 


Liv. Tel. 


r' 


1800 


I'.'l 






i' 


900 


I'.'O 


2006 Feb 20 


Liv. Tel. 


r' 


1800 


1'.'3 






i' 


900 


1'.'2 


2006 Feb 22 


Liv. Tel. 


r' 


1400 


a'9 


2006 Aug 22 


Liv. Tel. 


r' 


1400 


a'9 


2006 Sep 18 


Liv. Tel. 


r' 


1400 


a'9 


2006 Sep 30 


Liv. Tel. 


r' 


1400 


a'9 


2006 Oct 21 


Liv. Tel. 


r' 


1400 


a.'9 


2006 Nov 13 


Liv. Tel. 


r' 


1400 


a.'9 


2006 Nov 22 


Liv. Tel. 


r' 


1400 


a'9 


2006 Dec30 


LIRIS/WHT 


J, 


300 


I'.'i 






K, 


250 


I'.'i 


2007 Jan 8 


Liv. Tel. 


r' 


1400 


a.'9 


2007 Jan 14 


Liv. Tel. 


r' 


1400 


a.'9 


2007 Jan 21 


Liv. Tel. 


r' 


1400 


a.'9 



Table 1. Observing log 



due to atmospheric diff'erential refraction are very small, always 
below few percent, despite slit was not oriented at the parallactic 
angle. The flux calibrated spectra are presented in Fig.[T] 



2.2. Near Infrared photometry 

We obtained images in the J and filters using LIRIS at 
the same epochs as the near-IR spectroscopy. In both cases 
the images were obtained following a 5 point dithern pattern. 
Individual frames of 20 and 10 sec were taken in the J and Ks 
filters, respectively, giving as total exposure times those reported 
in Table|2l The images were reduced using the task Idedither in- 
cluded in the dedicated software package lirisdr, developed un- 
der IRAF by the LIRIS team. The main data reduction steps are 
sky subtraction, flat-field correction and finally a combination of 
the images after proper alignment. 

The photometric calibrati on was determine d by comparison 
with the 2-MASS catalogue dCutri et al.ll2003h . Zero points for 
our J and K, photometry were determined using about 20 stars 
with known 2-MASS magnitudes (brighter than / ~ 15.5) in- 
cluded in the LIRIS field. The typical dispersion between instru- 
mental and 2-MASS magnitudes is ~ 0.08. 
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Fig. 1. Final reduced spectrum of NRAO 150. The position of the most prominent features are marked. Deep atmospheric absorption bands are 
indicated by vertical bands. Top panel: Spectra corresponding to the bands Z+J observed in two epochs. The thicker red color line corresponds 
to 2005 Mar, and the thinner blue color line corresponds to 2007 Jan. Middle and bottom panels: Spectrum corresponding to the bands H and K, 
respectively. 
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Near IR 



Date 


JD 


J 


Ks 


1999 Oct 13 

2005 Mar 23 

2006 Die 30 


51465 
53453 
54100 


16.49±0.12* 
16.78±0.06 
16.52±0.04 


14.37+0.08* 
14.81+0.05 
14.54±0.05 


Optical 


Date JD 




V 





2005 Nov 2 


53677.45 


22.9±0.9 20.67+0.08 


19.43+0.10 


2005 Nov 11 


53683.4 


20.81+0.07 


19.72+0.04 


2006 Feb 3 


53770.35 


20.84+0.04 


19.19+0.02 


2006 Feb 17 


53784.35 


20.58+0.03 


19.32+0.05 


2006 Feb 20 


53787.48 


20.60+0.06 


19.68+0.05 


2006 Feb 22 


53789.38 


20.70+0.05 


19.57+0.06 


2006 Aug 22 


53970.65 


20.72+0.10 




2006 Sep 18 


53997.51 


20.88+0.02 




2006 Sep 30 


54009.56 


20.60+0.05 




2006 Oct 21 


54030.46 


20.62+0.01 




2006 Nov 13 


54053.39 


20.68+0.02 




2006 Nov 22 


54062.44 


20.78+0.02 




2007 Jan 8 


54109.38 


20.85+0.01 




2007 Jan 14 


54115.35 


20.67+0.04 




2007 Jan 21 


54122.40 


20.81+0.09 





Table 2. NRAO 150 - Optical near-IR Photometry. 
* Values corresponding to 2-MASS PSC. 
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Fig. 2. Identification charts of NRAO 150 in I and V fil- 
ters. The circle and square symbols mark NRAO 150 and the 
star 2MASS J03592889+505754, used as reference for spec- 
troscopy. Dashed square represent the plotted area in V band. 
North and East are toward the top and left of the frame, respec- 
tively. Numbers note the calibration stars proposed for this field. 



2.3. Optical Photometry 

In addition to the near-IR observations, we also obtained im- 
ages in the optical range using several telescopes at different 
epochs. The first observations were performed using the CCD 
camera mounted on the IAC-80 at the Teide Observatory during 
the nights of 2005 Nov 2 and 1 1 . We used filters V, R and I (see 
Table[TJ. Our target was detected in all bands (see Fig. |2]for ex- 
amples of I and V images). Given the faintness of the source, 
a large r diameter telesco pe, the 2.5 m Liverpool robotic tele- 
scope ( ISteele et al.ll2004l ) was used for the remaining observa- 
tions. Filters Johnson V and Sloan r' and i' were used, although 
most data were obtained in r' only (see Table[T]). 

All data were reduced in IRAF following the standard proce- 
dure, i.e., bias subtraction, flat-fielding, and image combination. 
The photometric calibration was based on images of Landolt 
fields obtained at photometric conditions. The photometry was 



obtained using the SExtractor program (iBertin & Arnoutslll996h 
and the AUTOMAG parameter. The resulting photometric mea- 
surements are listed in Table |2l We estimated that the detection 
of faint objects in our images is complete down to V = 23.1 
(24.0) mag, / = 21.7 (23.0) mag and /' = 20.4 (21.5) mag for 
S /N - 5 (3) within the observed fields. 

Based on our photometric calibration, we propose seven ref- 
erence stars for r' and /' filters which can be used for future mon- 
itoring. These stars are marked in Figure|2]and their correspond- 
ing magnitudes are listed in Table |3] 



3. Results 

3.1. Near I R spectrum 

A very prominent feature is observed in the spectrum of 
NRAO 150 in the H band at 1.65 jum (see Fig. [TJ- A less in- 
tense, but notable, feature is also seen in the J band at ~ 1 .23 //m 
. These are identified with the Ha and Hfi lines, redshifted by 
z - 1.517 + 0.002 and then providing the first redshift determi- 
nation of NRAO 150. The corresponding luminosity distance of 
the source is <iL = 11-2 x lO'^Mpc, for a //„ = 71 kms"' Mpc"', 
Q,„ = 0.27, and = 0.73 cosmologjQ. 

We measured the center, width and flux of the narrow and 
broad components of Ho- and Hp by fitting two Gaussian com- 
ponents plus a slope as the continuum (see Table |4]i. We started 
by fitting Ha since this part of the spectrum has higher signal-to- 
noise ratio and is not contaminated by other features such as Fell 
emission. The Ha profile is well fitted by a narrow component 
of FWHM=170A (1458 km/s rest-frame) and a broad compo- 
nent of FWHM=666 A (5745 km/s rest-frame). The broad com- 
ponent is blueshifted with respect to the naiTow component by 
62A (532 km/s). In the case of H/3 the S/N ratio is much lower 
and was not possible to let vary the line widths. Instead they 
were constrained from the best fitting value of Ha. 

Other broad features showing lower S/N are observed around 
4300 A (rest frame), corresponding to the Balmer limit, and at 
4450-4700, 4924, 5018 and 5150-5350 A, associated to blends 
of Fell emission (see Fig. [TJ. In order to remove and confirm the 
importance of Fell emis sion we have subt r acted the empirical 
Fell emission template of lBoroson & Greenl([T99l . scaled to the 
intensity of the observed features (see Fig[3]i. This template was 
generated from a spectrum of PG 0050-1-124 (IZw 1), which is 
well known by the strength of its Fell emission and t he narrow- 
ness of its H recombination lines (IOke&Lauel^ll979h . The tem- 
plate (kindly provided by Dr P. Marziani) was prepared by re- 
moving the lines which are not associated with Fe II transitions. 
It can be seen from Fig. |3]that after the Fell template subtraction 
the residual spectra show uniquely the broad H/3 feature, plus 
a shallow bump at 4300A corresponding to the Balmer limit. 
Given the radio loud AGN properties of NRAO 150, it is ex- 
pected to observe intense narrow forbidden emission lines such 
as [OIII] A 5007 A, [Nil] A 6548A, A 6584A, and [SII] AA 6716, 
6731 A. However, none of these features are detected at a sig- 
nificant level over the continuum. The implications of weak or 
absent typical NLR features will be discussed below in section 



* The cosmology calculator available in the Web 
|http : //mm . astro . ucla . edu/~wright/CosmoCalc . html| i 
was used. 
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ID 




K.A. 


Dec (J2UUU) 


r' 




/' 






(UJ:5y:ss) 


(+50. . ) 








1 


J03593205+5057513 


32.05 


57:51.4 


20.12 ± 


0.04 


19.10 ± 0.07 


z 




O J.O / 


J 1 .OL.L 


90 1 9 -t- 
ZU. i Z It 


u.uz 


1 S zLS -1-0 01 
1 O.'+O It U.U i 


3 




31.86 


57:45.0 


20.86 ± 


0.10 


19.79 + 0.14 


4 


J03592824+5058196 


28.24 


58:19.6 


19.62 ± 


0.01 


18.66 + 0.07 


5 




56.58 


57:47.3 


21.00 ± 


0.07 


19.75+0.14 


6 


J03592625+5057054 


26.26 


57:05.1 


18.98 ± 


0.01 


18.04 + 0.09 


7 


J035928 14+5058277 


28.18 


58:27.7 


18.76 + 


0.01 


17.70 + 0.04 



Table 3. Photometry of reference stars near NRAO 150. 
This table could be published only in the electronic version. 




4000 4200 4400 4600 4800 5000 5200 
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Fig. 3. Top panel: The two NRAO 150 spectra (histogram style line) 
taken at different epochs are represented together with a scaled version 
of the IZwl template (green continuous thin line). Bottom panel: The 
spectra after subtraction of the template. It can be seen that the resid- 
ual spectra consist basically on a broad H/J feature, a shallow bump at 
4300A corresponding to the Balmer jump, and very weak [0111] lines. 



3.2. Optical and near IR variability 

Our optical and near-IR photometric measurements span a range 
of almost 2 years (see Table |2]i, which is typically a sufficiently 
long time range to detect int r insic v ariability in radio loud AGN 
(e.g.. iTerasranta et al.ll200l l2005h . Indeed, NRAO 150 shows 
a nearly mon othonic radio flux inc rease by a factor ~ 4 from 
1991 to 2005 ( ITerasranta et al .12004 12005). However, we cannot 
make any firm claim about optical variability based on measure- 
ments in the r' band, by far the band with the longer and denser 
time coverage. The maximum magnitude variations observed are 
~ 0.2 mag, which are not larger than 3 times the typical measure- 
ment uncertainty. It is worth to note that part of the data were 
taken using different combinations of filters, cameras, and tele- 
scopes which does not guarantee their homogeneity. Also, the 
faintness of the source in the optical range prevents to obtain 
more accurate photometric measurements. 

As in the optical range, our near-IR photometric measure- 
ments (performed in 2005 March and 2006 December) do not 
show significant variability with amplitude above ~ 0.3 mag. 
(~ 0.25 or 30%). A similar result is obtained when our near-IR 
measurements are compared with those of the 2-MASS survey 
in 1999 (see Table |2]i. In this case, maximum variability ampli- 
tudes of AJ = 0.3 mag and AKg, = 0.4 mag are found, which 
contrasts to the factor of 2 flux density increase reported by 
ITerasranta et al.l ( l2005l) at radio wavelengths for the same time 
range. This suggests that the process responsible for the radio 



emission is not connected to the one responsible for the optical 
and near IR emission. 

4. Discussion 

4.1. Black hole mass and its accretion efficiency 

We use the empirical relationships provided by 
IVestergaard & Peterso nl (l2006l) to estimate the central black hole 
mass in NRAO 150. The black hole (BH) mass can be estimated 
from the luminosity and width of the broad componen t of th e 
Hy6 line using expresion 6 in IVest ergaard & PetersonI ( |2006|) . 
Given the low S/N ratio of the H/? line detected in our spectral 
measurements, the use of the line luminosity is preferred against 
the value of the co ntinuum at 5100A (used in expression 5 of 
IVestergaard & Pete rson 2006). The measured broad component 
luminosity is 5.6 ± 1.2 x lO^^ergs^^ which increases up to 
2.3 + 0.5 X lO'^'^er^ i"' after applying the galactic extinction cor- 
rection derived in this work (see Appendix|A]i. Combining these 
values with those for the line width (see Table |4]i the resulting 
black hole mass is 1 .95 X lO'^Mo, which converts to 4.68 x IO'^Mq 
after extinction correction. Using the estimate of the black hole 
mass we can also determine the corresponding Eddington 
luminosity [LEdd = 3.3 x IQ'^{Mbh/Mq) Lq], which results in a 
value of 1.54 x IO'^^Lq. In order to compute the accretion rate, 
this value must be compared with the bolometric luminosity 
of the accretion disk, which results in Ldisk - 7.8 x 10 '^Lq, 
and Ldisk - 46.9 x lO'^L© after extinction correction. The disk 
luminosity was computed by integrating the measured fluxes 
within the optical-UV range in the rest frame, this value is 
multiplied by 27rD- to obtain the luminosity. Here it is assumed 
that the radio-mm and the X-ray emission are related to the 
relativistic jet. This is justified given the uncoupled variability 
of the optical-UV spectral ranges with regard to the radio-mm 
one as reported in Section 13.2! Thus the Eddington luminosity 
ratio is Ldisk/LEdd = 0.30. 

These values can be compared to those previously reported 
in the literature for similar objects. The black hole mass of 
NRAO 150 is above the hi ghest masses found in the li tera- 
ture for low redshift AGNs (IVestergaard & Petersonll2006l). al- 
thoug h the estimated value for M 87, is 3.4 X lO^M^ (iGrahamI 
l2007l) . comparable to the one of NRAO 150. In contrast, flie 
black hole mass of NRAO 150 is well in the range of typical 
masses of luminous high redshift AGNs (She r pmer et al. 20()^ . 
Several works (e.g., iMcLure & Dunlopll2004l: rShen et al.il2008l) 
have tried to measure the dependence on redshift of the mean 
BH mass. They reach a common result: the mean BH mass in- 
creases with redshift, although this dependence is dominated by 
the M almquist bias (Vestergaard 2008j). Recently, Labita e t all 
(l2009l) claimed that the maximum BH mass evolves with z as 
log{MBH{max)/MQ) ~ 0.3 z -1-9, whereas the maximum Eddington 
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Center 


FWHM 


Flux 


Lum. 


Line 




(A) 


(A) 


km/s 


Obs. 


Ext. Cor. 


Obs. 


Ext. Cor. 


Ha(n) - 


2005 March 


16501. 1±0.7 


72±1.7 


1458 


22.0±0.5 


51.7+1.2 


32.9+0.7 


77.1 + 1.8 


HQ-(b) - 


2005 March 


16439±7 


283±16 


5745 


21.7±0.7 


51.0+1.6 


32.4+1.0 


76.1+2.4 


Hj8(n) - 


2005 March 


12255±5 


54±1.2 


1458 


2.97±0.3 


12.0±1.2 


4.4+0.4 


17.9+0.6 


mb)- 


2005 March 


12209±5 


210±12 


5745 


3.74±0.8 


15.1+3.2 


5.6+1.2 


23+5 


H/3(n) - 


2007 Jan 


12254±9 


54±1.2 


1458 


2.14±0.4 


8.6+1.6 


3.2+0.6 


13+2 


H8(b)- 


2007 Jan 


12208±9 


210±12 


5745 


3.45±0.9 


13.9+3.6 


5.1+1.3 


21+5 



Table 4. Flux units are 10 '^ergs 'cm ^. Luminosity units are lO'^^'ergs 



ratio ( 0.45) is essentially constant with z. However. iLabita et alJ 
(|2009|) suggest that these results are unaffected by the Malmquist 
bias. Our estimate for the BH mass in NRAO 150 is slightly 
above the maximum value, and the Eddington ratio is also close 
to the maximum value, which indicates that NRAO 150 pos- 
sesses a very massive and efficiently accreting BH. 

4.2. Fell emission and weakness of the NLR emission 
features 

As reported above, the spectrum of NRAO 150 presents promi- 
nent H recombination lines plus intense Fell emission, in 
contrast with very weak or absent lines charac t eristic of the 
NLR, such as [OIII] (see Fig. O. iNetzer et al.l (|2004 found 
that about one third of very high luminosity AGNs do not 
show strong [OIII] lines. In addition, an anticorrelation be- 
tween EW(FeII)/EW(Hjg) and EW(rOIIII) (iBoroson & Greed 
[T99I r^ian"& Willsl I2OOI INetzer et all |2004 is reported. In 
our NRAO 150 spectra, the estimated EW for H/3, Fell and 
[OIII] are around 55, 60 and 6, although the uncertainties are 
large given the difficulty to deblend the spectral features. These 
values are consistent with the mentioned results. Using these 
values we checked that NRAO 150 is placed in Fig. 9 of 
INetzer et al.l (|2004|) at the locus of high luminosity and high- 
re dshift QSO sample The location of NRAO 150 in Fig. 3 
of INetzer etaP (12004 is intermediate between the one of PG 
OSOs ( lBennert etal.ll2002h and z > 2 QSOs, according to its 
Hy6 luminosity. Netzer et aP (l2004 suggest that for highly lu- 
minous QSOs the NLR becomes extremely large (even larger 
than the size of the largest galaxies) as expected from the "nat- 

1 /2 

ural" Rnlr i^ig,, dependence, which makes NLR to disap- 
pear becoming dinamically unbound. The equivalent size of 
the NLR in NRAO 150, as predicted from the HIS luniinosity 
(~ 4 X lO'^^ergs-i), is about 20 kpc (Fig. 3 in lNetzer et al]|2004 . 
which is well above the sizes commonly reported for typical 
NLR, hence explaining the weakness of the NLR emission fea- 
tures in our spectra. 



4.3. The spectral energy distribution 

Fig. ID shows the spectral energy distribution (SED) of 
NRAO 150. In addition to the photometric data presented here, 
we searched for data covering the widest possible wavelength 
range, from radio to X-rays. At radio bands, the lowest frequency 
measurements, obtained at 4, 8.4, 22, and 43 GHz with the VLA 
on March 2005, were taken from the National Radio Astronomy 
Observator}0 data archive. The next three measurements were 
observed at 86, 142, and 229 GHz under the general IRAM 30- 



meter Telescope AGN Monitoring Program dReuter et al.ll 19971 
and references therein) on March 2005 (August 2005 for the 
larger frequency one). The optical-UV rest frame observations 
corresponds to the near-IR and the optical observations pre- 
sented here, and were acquired on March 2005 and November 
2005, respectively. On this spectral region, empty triangles in- 
dicate the extinction uncorrected optical and near-IR measure- 
ments. Empty circles symbolize such measurements corrected 
for the Ga lactic extinction estim ated from IRAS far-IR emis- 
sion maps (ISchlegel et al.lll998l) . whereas filled circles indicate 
the same measurements corrected for the more accurate Galactic 
extinction estimated here (see Appendix |A]i. The X-ray data, 
acquired by ROSAT from August 1990 and Fe bruary 1991, 
were c orrected by Galactic extinction as reported in lAgudo et al.l 
(l2007l) . For the plot in Fig. H a photon index F = 1.7 + 0.1 
was assumed based on the synchrotron spectrum measured from 
the IRAM 30 m observations presented here at 86 and 142 GHz 
(spectral index a = 0.7 + 0.1). 

The SED of NRAO 150 shows two prominent bumps: 
one peaking around 1 mm and another one at 2000 A. The 
low frequency bump observed on the SED resembles those 
typical of high power flat-spec t rum radio-loud sources (e.g ., 
iGhisellini. Tavecchio. Ghirlandal l2009t iGhiseUini etaP l2009l) . 
This first bump is attributed to strong Doppler-bosted syn- 
chrotron relativistic -jet-emission. Also the X-ray domain of 
the SED is typical from inverse Compton emission from the 
jet, whereas the large luminosity optical-UV bump (L ^ 
10''^erg/s at V » lO'^^'^Hz) is not so commonly reported 
for this kind of sources. A prominent peak at optical-UV 
(rest frame) wavelengths is typical of Seyfert galaxies and 
is thought to be produced by thermal emission from the 
AGN accretion disk. However, there is an increasing num- 
ber of high-power radio-loud blazars for which this emission 
feature is being reported (e.g, iRaiteri et al. 2007 : AbdoetalJ 



20091: iD'Ammando et al] 120091: iGhisellini & Tavecchio 



I2009t 



Gh isellini. Tavecchio. Ghirlandal 120091: IGhisellini et al.l I2OO9I1 
Indeed, nearly h alf of the radio loud blazars in the sam- 
ple considered by IGhisellini. Tavecchio. Ghirlandal (I2009h were 
recently reported to show a similar optical-UV large ex- 
cess. There is a consensus to consider this emission to 
come from thermal emission from the accretion disk (e.g , 
IRa iteri et al.l 120071: lAbdo et all 120091: Ib'Ammando et al.l 120091: 



^ The National Radio Astronomy Observatory is a facility of the 
National Science Foundation operated under cooperative agreement by 
Associated Universities, Inc. 



Ghi sellini & Tavecchio I 2009t IGhisellini. Tavecchio. Ghirlandal 
i2009l: IGhisellini et al.ll2009l) . that is not outshined by the syn- 
chrotron bump from the jet when such bump peaks at v << 
10'^^ Hz. We also agree on such expl anation for the case o f 
NRAO 150 for the reasons outlined bv IGhisellini et al.1 (l2009l) . 
Moreover, this also explains the above-reported lack of variabil- 
ity in the observed optical and near-IR bands of NRAO 150, 
whereas the radio and mm spectral ranges showed a factor 5 long 
term variability. This is easily explained if such observed optical 
and near-IR (rest frame optical-UV) emission is produced in the 
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Fig. 4. Spectral energy distribution (SED) of NRAO 150 from the ra- 
dio to the X-ray spectral range after correction for galactic extinction. 
Filled circles represent Galactic extinction estimated in this work (see 
Appendix [At. instead empty circles represent ext inction estimated from 
IRAS far-IR emission maps dSchlegel et aljl993) . Empty triangles indi- 
cate the extinction uncorrected optical and near-IR measurements, note 
the dramatic change of the SED after extinction correction is applied. 
The measurements presented here were obtained essentially contempo- 
raneously during year 2005 (see text), except for the ROSAT X-ray data 
(bow tie), that were acquired between August 1990 and February 1991. 



accretion disk and not in the relativistic jets where the extreme 
variability comes from. 

4.4. Source classification 

From the VLBI point of view, NRAO 150 shows a one-sided 
powerful relativistic-superluminal jet typical of blazars, which 
indicates that the axis must be aligned close to the line of sight. 
Moreover, its radio spectrum is typically flat, which classifies 
the source as a flat spectrum radio quasar (FSRQ). The opti- 
cal rest frame spectrum shows intense and broad (~ 5000 km/s) 
Hq. and emission, so the source is a type 1 AGN, consistent 
with the disk axis being aligned close to the li ne of sight accord- 
ing to the standard AGN Unification scenario dUrrv & Padovanil 
Il995h . The optical and radio properties are both consistent rel- 
ative to the orientation of the radio jet and the accretion disk. 
The fact that the forbidden [OIII] lines are weak points to a very 
luminous central engine, which is also consistent with the high 
Eddington rate determined from the optical/UV luminosity rela- 
tive to the BH mass. The large accretion rate is consistent with 
the large radio loudness of NRAO 150. The BH mass is also 
among the highest derived for quasars d^sterg aard & PetersonI 
I2006t I VestergaardI 120081: lOhiseUini et al.l l2009l and references 
therein). 



5. Summary and conclusions 

We have determined, for the first time, the cosmological distance 
of NRAO 150; one of the brightest radio to mm AGN sources in 
the northern sky; by means of its spectroscopic redshift (z= 1 .5 17 
or c/l - 11164Mpc). Given the low Galactic latitude of the 
source, its optical spectral lines have remained hidden to us for 
decades. The new near-IR spectra presented here have revealed 
very intense Ha emission together with less prominent HyS and 
Fell blended lines. 



The line width of Ho- points out that NRAO 150 is a broad 
line AGN, whose disk axis must be oriented close to the line 
of sight, consistent with having a powerful one-sided superlumi- 
nal radio-mm jet. Based also on the radio to mm spectral and 
variability properties (including those observed with VLBI), we 
classify NRAO 150 as a FSRQ blazar. Consistent with previous 
observations of high redshift and highly luminous quasars, the 
optical rest frame spectrum of the source show weak or absent 
spectral features typical of the NLR -such as [OIII]/15007 A. 
This supports the idea that the highest luminosity accretion disks 
(and probably the most massive BH) favour the unbounding of 
their NLR by accretion disk radiation. 

Using empirical relationships between HfS line width and lu- 
minosity we estimate that the central engine in NRAO 150 pos- 
sesses a large black hole with mass ~ 5 x IO^Mq. The radio 
to X-rays SED of the source shows two prominent bumps: one 
peaking at millimetre wavelengths -typical of synchrotron radi- 
ation from high power blazar jets-, and another in the near-UV 
(~2000A) that is attributed to thermal emission from the accre- 
tion disk. The good spectral coverage of the disk emission allows 
us to make a reliable measurement of the bolometric luminos- 
ity of the disk, which turns out to be accreted at ~ 30 % of the 
Eddington rate. 

This, together with the large BH mass of the source, its 
prominent BLR line luminosity, and its highly luminous syn- 
chrotron spectrum sets NRAO 150 among the most pow- 
erful FSRQ blazars. Such sources are also the most lumi- 
nous hard X-ray and y-ray inverse-Compton blazar emitters 
(Ghisellini & Tavecchio 2009; Ghisellini, Tavecchio, Ghirland^ 
I2OO9I: iGhisellini et al.l |2009), and are routinely monitored by 
high energy space observatories as Fermi. We predict that 
NRAO 150 is one of such sources. However, it has not still been 
detected in y-rays. The low Galactic latitude of the source im- 
plies a challenge for y-rays observatories to detect it. If that is 
possible in the future, modeling of the whole broad-band spec- 
trum SED would allow for further investigation of the intrinsic 
physical parameters of this powerful blazar 
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Appendix A: Galactic Extinction correction 

We have already mentioned that the low galactic latitude of 
NRAO 150 and its associated high galactic extinction has de- 
fied the optical identification and studies of this radio-source. In 
order to compute the intrinsic properties of NRAO 150 we have 
to correct for the galactic extinction. At a first approach, we have 
taken the extinction from SIMBAD database, which is estimated 
from far infrared em ission maps build by combining IRAS and 
COBE/DIRBE data dSchlegel et alJll998h . The quoted value is 
E(B-V)= 1.474, which implies Ay = 4.5 to Ak = 0.54. 

In order to have another estimate of the galactic extinction 
towards the direction of NRAO 150, we built a color-magnitude 
diagram (Fig. lA.ll ) using the near infrared colors (J and K from 
2MASS database Cutri et al. 200l) of neighbouring stars. The 
color magnitude {J-KsJ) diagram stellar distribution was com- 
pared wothjhetheoreticalisochrones for 2MASS filters retrieved 
from iGirardi et alj (|2002|) . We intend to recognize part of the 
Main Sequence in the color-magnitude diagrams. Indeed, we 
recognized two possible Main Sequences (see Fig. lA.ll ). one of 
them, more crowded, showing an excess J - K ~ 0.7 and a sec- 
ond sequence, less populated, with J - K ~ 1.3. We identified 
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Fig. A.l. Color-magnitude diagram of a field of 30 x 30 arcmins, 
centered on NRAO 150. The two lines correspond to the theoretical 
Main Sequences of an isochrone with logf = 9.1 Gyr and z=0.019 of 
iGirardi et al.l l l2002h . Continous line fits the more distant and extincted 
stellar population while dashed line represents the fit to a less extincted 
population of stars. 

the first population to be in front of the Galactic disk and the 
other one to be behind the disk and presenting a higher galactic 
extinction. 

If we consider the solar metalicity, z=0.019, we found that 
the best possible fits for both stellar populations correspond to 
isochrones with logf - 9.1 Gyr (see figure IATTT i. For these fits, 
we compute E(J-K)=I.O + 0.1 for a stellar population located at 
distance module of 10.2, and E(J-K)= 0.12 + 0.3 and a distance 
module of 9.5 for the less extincted and closer population. We 
are only interested in the largest value of E(J-K), which gives a 
value E(B-V)=1.7 + 0.2 (E(B-V)=1.69 E(J-K)). This value is in 
reasonable agreement with the estimate provided by SIMBAD 
database. 



